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Abstract. A family of non-equilibrium statistical operators is introduced which differ by the system age 
distribution over which the quasi-equilibrium (relevant) distribution is averaged. To describe the nonequi- 
librium states of a system we introduce a new thermodynamic parameter - the lifetime of a system. 
Superstatistics, introduced in works of Beck and Cohen [Physica A 322, (2003), 267] as fluctuating quan- 
tities of intensive thermodynamical parameters, are obtained from the statistical distribution of lifetime 
(random time to the system degeneracy) considered as a thermodynamical parameter. It is suggested to 
set the mixing distribution of the fluctuating parameter in the superstatistics theory in the form of the 
piecewise continuous functions. The distribution of lifetime in such systems has different form on the dif- 
ferent stages of evolution of the system. The account of the past stages of the evolution of a system can 
have a substantial impact on the non-equilibrium behaviour of the system in a present time moment. 
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1 Introduction 

One of the most promising ways of development of the 
description of the non-equilibrium phenomena is provided 
by the method of the non-equilibrium statistical opera- 
tor (NSO) P^[5]. In [5 a new interpretation of the NSO 
method is given, where the operation of taking the in- 
variant part [IJ or the use of the auxiliary "weight func- 



tion" (in the terminology of [MS]) in NSO are treated 
as averaging of the quasi-equilibrium statistical operator 
over the distribution of the past lifetime span (age) of 
the system. In [5J|5] it is noted, that multiple choice of 
the "weight functions" can be taken. In [5j a uniform dis- 
tribution over an initial moment Iq is considered, which 
after the change of integration order reduces to the expo- 
nentially distributed weight function Pq{u) — exp{— ew} 
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in Such distribution is the limiting case of the life- In [2 a physical interpretation of the exponential dis- 

time distribution 7J, that is the distribution of the first tribution for Pq{u) is given: the system freely evolves as 

passage time of a given level. The term "lifetime" denot- an isolated system governed by the operator of Liouville. 

ing the time of the first passage of a level was used in In addition, the system undergoes random transitions at 

[7]. Encountered in the literature is also the term "non- which its phase point representing the system spreads 

equilibrium relaxation time" , as well as some others. from one phase trajectory to another one in a random 

The form of the function pq{u) = exp{-eu} in ^ fashon with an exponential probability under the influence 

is connected with the form of the source in the Liou- of a " thermostat" , thus average intervals between random 

viUe equation for NSO. In ^ the sources in the Liou- J^mps increase infinitely. This feature reflects in the pa- 

viUe equation different from that introduced in the NSO ramctcr of the exponential distribution tending to zero 

method in ^ are considered. The approach of the present after the thermodynamic hmiting transition. Real physi- 

paper differs from the methods used in But the use ^^1 systems have finite sizes. The exponential distribution 

of the distribution of the system lifetime in the present describes completely random systems. The influence of the 

work can be compared with that in ^ enlarging the set of surrounding on a system can have organised character as 

macroobservables; besides the common physical macroob- ^^11, for example, for the systems in the non-equilibrium 

servables this approach includes additional ones, namely steady-states with input and output flows. The ways of 

the life span (lifetime) ^^"^ interaction of the system with surrounding can be dif- 

__ ferent, therefore various choices of the functions p„(u) are 
In [lOj other alternative derivations of NbO are per- 

justified. 

formed, following the ideas proposed by McLennan [IT], 

and a relation with an earlier proposal by I. Prigogine [12j 

is discussed. The source in the Liouville equation can in Nonextensive statistical mechanics [HHS] can be re- 
principle take different forms. The form of a source used in ^^^^^^ embedding of the common statistical me- 
[TUH is a specific case which can be obtained under the chanics into a more general framework. Many complex 
assumption of the weak coupling limit of the interaction gygtems exhibit a spatio-temporally inhomogeneous dy- 
of the system with its environment. namics that can be effectively described by a superposition 
In [13j it is shown, what is the impact of changing the of several statistics on different time scales, termed "su- 
system lifetime distribution to the non-equilibrium prop- perstatistics" |16[|17) . Nonequilibrium situations are de- 
erties of system with finite volume. In the present work we scribed by a fiuctuating parameter (3, which can be, for 
consider infinitely big systems with infinite average life- example, the inverse temperature. The generalization of 
times as well. the Boltzmann factor exp{— was introduced in the 
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following form: 2 Modifications to the nonequilibrium 

description 

POO 

B{E) = dfi' f{l3') exp{-/3'£;}. (1) a new interpretation of the method of the NSO is 

given. Let us consider now, what consequences follow from 
such interpretation of NSO. Setting various distributions 
The type of superstatistics induced depends on the for past lifetime of the system ^^(m), we obtain a family of 
probability distribution /(/3) of the parameter (}. The spe- NSO, where the exponential distribution in Zubarev NSO 
cial case of these superstatistics, with the function / in the is a particular choice of the form of the (arbitrary) lifetime 
form of gamma-distribution, appears in the nonextensive distribution Pq{u). The class of NSO from this family can 
statistical mechanics |16[|171I141[T5] , describing a number of be related to the class of the lifetime (or age) distributions 
physical phenomena which are not satisfactory described (taken, for example, from the stochastic theory of storage 
by the Boltzmann-Gibbs statistics. In the present work the processes, the theory of queues etc) and to the relaxation 
superstatistics like (P) (together with its generalization) is properties of that class of physical systems which is inves- 
obtained starting from the distribution which contains a tigated. The general expression for NSO with an arbitrary 
lifetime of a statistical system as a thermodynamic vari- distribution pq {u) is 
able [I8 l ll9 p 20 p 21j. Section 5. This distribution has been 

hiQ{t)— I pq{u)\n.Qq{t ~ u,—u)Au~ (2) 
applied earlier to the description of aerosol behaviour [22] , " ^ 

I n\ f f f f \ A \ ~ "^-^ A 

and neutrons in a nuclear reactor [23]. = in ^^(t, 0) - / I / Pq[u) du j — du , 



In gq{t,t,) = e{-*i^/*^hne,(t,0)e{*i^/*^}; 
<?(t) =ln Spexp{-^F„(i)P„}, 



In several works, for example in [TBlflT] . the distribu- \ngq(t,0) = -^(t) - ^F„(t)P„; 

n 

tion function f{(3) is introduced as some continuous func- 
tion expressed through arbitrary analytical form of the 
distribution of a random variable /3. But the definition of 

the continuous density of distribution assumes its piece- where H is Hamiltonian, In g(t) is the logarithm of the 

wise continuous character when the density of distribution NSO, In gq{t, 0) is the logarithm of the quasi-equilibrium 

has finite number of breaks. Real nonequilibrium systems, (or relevant) distribution; the first time argument indi- 

as a rule, are spatially non-uniform. This behaviour can cates the time dependence of the values of the thermo- 

be mathematically described by the piecewise continuous dynamic parameters Fm', the second time argument t2 in 

functions, the examples of which are given in the present ^?g(^l, denotes the time dependence through the Heizen- 

work. berg representation for dynamical variables Pm on which 
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gq{t,0) can depend [IlEllSlllllG] , integration by parts in role of the form of the source term in the Liouville equation 
time is carried out at J Pq{y) dy\y=n ~ —1; J Pqiu) dy|y_+oo = in^ NSO method has never been investigated. 
0. If Pq{y) = £exp{— ey};£ = l/(r), then the expression Let us notice, that in the case when the value 
([2]) passes in NSO from [1,2 . In [4] the auxiliary weight dlTigq{t — y, ~y)/ dy (the operator of entropy production 
function [SJIS] Pq{u) — e exp{— ew} was interpreted as the a [Ull]) in the second term of the r.h.s. ([2]) docs not depend 
density of the probabihty distribution of the lifetime F oi on y and can be taken out from the integration on y, this 
a system. There -T is a random variables of the lifetime second term takes on the form <j{r), and the expression 
of a system from the moment to of its birth till the cur- thus does not depend on the form of the function Pq{y)- 
rent moment t; e^^ = {t ^ ^o); {t — to) — (F); u = t — to- It is the case, for example, if Qq(t) ^ exp{— crt}, a — const. 
This value represents the age of system. The operator of In |25| such a distribution is obtained from the principle of 
internal time describing the age of a system was also in- maximum of entropy applied to the set of average values 
troduced in I.Prigogine's works (for example, see [T^). If of fluxes. 

the interval t~to — uis large enough (that is the details of As it is known (for example, [7]), the exponential dis- 
an initial condition as dependence on the initial moment tribution for lifetime 

to are nonsignificant and nonphysical [Tl[5] ) , it is possible , , , , , . , 

^9(2/) = eexp{-e2/}, (4) 

to introduce the minimal lifetime F,-ain — A and to in- 
used in the works of Zubarev [HISj is the limiting distri- 

tegrate in ^ over the interval (1^1,00). It is possible to 

bution for the lifetime, valid for large times. Thus, in [1] 

specify many concrete expressions for lifetime distribution 

[2] the thermodynamic results are obtained, which in this 

of a system, each of which possesses its own advantages. 

context are valid for all systems. 

Each of these expressions induces some form of a source 

For the NSO ^ with the function (g]) in the form 

in the Liouville equation for the nonequilibrium statisti- 

suggested by Zubarev the value in the second item is 

cal operator. In the general case of an arbitrary function 

Pq{u) the source is: " j Piiv) ^2/ = exp{-£2/} ^ I ~ ey + [eyf /2 - ... = 

J ^Pq{Q)h^gq{t,Q)+ r ^PM(^inQq{t-y,^y))dy 

' (3) l-y/{F)+yy2{Ff-.... 

(if the value Pq{0) diverges, it is necessary to choose the Evidently the average lifetime tends to infinity, (F) 00, 

lower limit of integration equal to some Fmin > 0). Such and the correlation ([5]) tends to unity, 
approach corresponds to the dynamic Liouville equation Besides the exponential density of probability ([3]), the 

in the form of Boltzmann-Bogoliubov-Prigogine [51151112] , Erlang distributions (special or general form) , gamma dis- 

containing dissipative items. In [24] it was noted that the tributions etc (see [^51l27| ). as well as the modifications 
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considering subsequent composed asymptotic of the de- tems with infinite lifetime, for example, for the systems of 

composition 28J can be used as candidates for the den- infinite volume after thermodynamic limiting transition. 

sity of lifetime distribution. Gamma distributions describe » j ± ± J.^ i. • i- jri ■ 

Amendments to the unity term m the equation (|5|) m 

the systems whose evolution acquires some stages (num- r™ , • u- i u ^, j-u ■ r-n j. 

^ & V become vanishmgly small when the size ot the system 

ber of these stages is given by the order of gamma dis- , •, rr j 4- • c v • -ti, j 1 

° ° " and its average lifetime tend to mhnity, as m the model 

tribution). Considering actual stages in non-equilibrium j- . -t m j ■ r? 1 ato^ ^ j. r 

' t> t> ^ distribution ([4]) used m Zubarev NbU. hor the systems ol 

systems (chaotic, kinetic, hydrodynamic, diffusive and so n •, • 1 , ,-n i j- ^ 

■' ^ T J J hnitc size but still exponential distribution these terms re- 

forth 1291 1. it is possible to justify the use of gamma dis- • j j_ ^ • jfi rm 

' — J J o g^l^ jj-^ nonzero amendments to the expression (|5|). ihus, 

tributions of a kind these additional terms to NSO and hence to the kinetic 

p (y) = £{ey)''~^ exp{—£y} / r{k) (6) equations, kinetic coefficients and other non-equilibrium 



{r(k) is gamma function, passing at fc = 1 to exponen- 
tial one H])), and to understand their importance in the 
description of non-equilibrium properties. The piecewise- 
continuous distributions corresponding to different stages 
of evolution of the system will be used below. 

More accurate specifying the shape of the function 
Pq{u) in comparison with limiting exponential function 
([U allows to describe in more detail the real stages of 
evolution of a system. Every form of the lifetime distribu- 
tions has certain physical sense. In |13| additional terms 
to NSO in Zubarev form for the gamma distribution ^ 
are obtained. 



characteristics of the system are in fact the impact of the 
finitcness of size and lifetime of the system. They do not 
result merely from the choice of particular form of the 
lifetime distribution in the system. In what follows we 
investigate whether there exist such distributions of the 
lifetime of system for which the additional contribution to 
NSO differs from Zubarev NSO even for infinitely large 
systems with infinitely large lifetime. 



We will consider several examples of choosing the func- 
tion Pq{u) in ([2]). We limit ourselves to the piecewise- 
continuous distributions, from where the results different 
from [1,2J follow. There are numerous experimental evi- 
dences of such changes of the distribution of lifetime of 
the system pq (u) on the time scale of the life span of the 
system. In [501151] the transition of the distribution of the 
The amendments to NSO in the Zubarev form in the Sec- first passage processes from Gaussian regime to the non- 



3 Systems with infinite lifetime 



tion 2 and in [13j are obtained for the systems of finite Levy behaviour in a specific time moment is shown. Real 

sizes and lifetimes. We will show further, how the same systems possess finite sizes and finite lifetime which im- 

effects, involving the infiuence of the past of system on plies the influence of surrounding on them. The fact that 

its present non-equilibrium state have impact on the sys- these sources do not vanish in a limit of infinitely large 



V.V. Ryazanov: Lifetime distributions in the methods of non-equihbrium statistical operator and superstatistics 



systems, is related to the openness of the system, hence 
to the influence of its surrounding. 
3a) We shaU set 

ka'' 



u < c: 



(7) 



6eexp{— eu}, u>c. 



From the condition of the normalization Pq{u) = 1 
the normalization multiplier is 6 = e^'^ \ — ■ — ) . Average 



a + c 



lifetime (F) for the distribution ([7]) is equal to 



in- 



k - 1 



a 



a + c 



1 , ^ (kc + a) 



(8) 



The value ^ (r) ^ oo at e ^ 0. From ^ we find 

fc 



Pqiu) 



u < c: 



+ a ^ 

6exp{— eit} , u > c. 



(9) 



The source in the right part of Liouville equation for 
the distribution ([7]) in accordance with the expression ^ 
equals 

fc(fc + l)a'= 



k 

J = -ln£-g(t,0)- / 
a Jo 



{u + a) 



k+2 



S du- 



a + c 



e^e-""5du, 



where S = \ngq{t - u, -u). The distribution of NSO © 
in the case of ([7]) is 

ln£-(t) = \n,ub g{t) + A; 



a 



a + u 



a + c 



adu+ 
e~^"crdM, 



where a = din g{t — u, —u) / du; hizub Q{t) = higq{t,0) + 
Jq°° e~'^"cr du is the distribution obtained by Zubarev in 
[1112], and Z\ is a finite amendment to it. 



3b) We will consider now the distribution of the foL 
lowing kind: 

d , u < c; 

e exp{— £u} , u > c. 



Pqiu) 



(10) 



From the condition of the normalization we find 

d = - {1 — e^^'^). The average lifetime is 
c 

(r) = — + -e-"(l + ec). 



(11) 



The average lifetime {F) — > oo at e ^ 0. The source ^ 
in the Liouville equation in this case equals 



J = din Qq{t, 0) 



e e 



The amendment to the Zubarev form of NSO is 



A = 



'-n-e- 



-ec 



■ du. 



We see in this case, that an additional memory term 
of the system stems from the finiteness of its size, and the 
limited memory effect is observed. It is possible to consider 
other examples of the functions Pq{u) which result in the 
limited memory effects. 

3c) For the exponential density of distribution but 
with different intensities in different time intervals 

{eiexpj— Eim), u < c\ 
(12) 
&£2exp{-e2M} , M > c, 

from the condition of the normalization it follows that 
(r) = - [1 - e-^i=(l + £ic)] + -e-^^^(l + 820) . (13) 

£l £2 
At £2 0, {F) oo. 

J = ei\Tigg{t,0) - [ £^e~^i"5du- 

"'0 



sle-'^^Sdu; 



V.V. Ryazanov: Lifetime distributions in the methods of non-equihbrium statistical operator and superstatistics 7 

electric field, as linear reaction on a mechanical pertur- 
c bation, that is we regard the electric conductivity in the 

/"OO 

/ [e^"^^" — l] crdu. linear approximation, following the results of [3 and, as 

Jo 

in [3], we limit ourselves to the important special case - 
the reaction of the equilibrium system to the spatially ho- 
mogeneous variable field 



The natural question is now why do the examples of 
this section differ from the examples of Section 2. In the 
interpretation of [2 it is the random value to in u ~ t ~to 

that fluctuates. In [2] the limiting transition is performed E'^{t) = f —e (tj). 

for the parameter e,e — > in the exponential distribu- 
tion Pq{u) = e exp{—eu} after passing to the thermody- 
namic limit. In the interpretation of this corresponds 
to the fact that the mean lifetime of the system (F) = where P is the operator corresponding to the vector of po- 



The Hamiltonian of perturbation is given by 
Hi = -PE°{t), 



(t — to) — 1/e oo. But average intervals between ran- larization of the system. In the coordinate representation 
dom jumps infinitely incease, exceeding the lifetime of the this operator is written as 



^ e^r, 



where Ci is the charge of the particle, and is its position 
vector. The operator of current is 



system. Therefore a source term in the Liouville equa- 
tion tends to zero. If the change of the distribution Pq{u) 
caused by the influence of the surrounding, occurs on the 
time interval of the life span of the system, as in examples 

3a)-3c), its impact remains even if the mean lifetime tends J = P = e f, = — n,- 

j j 

to infinity. 

where pj is particle momentum, m is mass of a particle. 

We choose a model in which the Coulomb interaction is 
4 Application of the distributions of Section 3 taken into account as a self-consistent screening of the 
to the conductivity field, i.e. we take E = E'^. The most essential difference 

from ^ consists in the replacement of the Laplace trans- 



formation used in [3' , that is 

/■oo 

{A; 5}^+,, = / dt e^("+-)*(A(i), B)), (e > 0) , (14) 



On the example of conductivity we will investigate, what 
are consequences of the change of the type of functions 
Pq{u) and g{t) as compared to the exponential law for 

Pq{u), used in [3]. where {A{t),B{t')) = dxTi-{AA{t)AB{t' + if3hx)geq} is 

The determination of the conductivity coefficient by the time correlation function [3], by the another integral 

the NSO method is considered in [321[351[Mll35j . In this transformation. So, for the example 3a) with the distri- 

section we will describe the transport of charges in the bution of the form ^ the expression (fT4|l is replaced 
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by For the distribution (fTU]) in 3b) the Laplace transform 

/ Q \ of the kind ([14]) is replaced by 

{A; B)^,a,k + e"" —— {A; B)^+,e;(c,oo) , (15) 
\a + cj 

-~d{A; B)^.c,t + {A; i?}cj+ie:(c,oo) , 

where 

k where {A;B)^.c,t = dt e"^H{A{t) , B) , the value d is 



{A-B)^,a,k^ f^dte-'{-^^ {A{t),B)- 



given in HU]), {A; B)^^^^.^^^^) is given in In- 

(A;B)^+i,.(e,oo) = / dte^'^+''^\A{t),B). (16) stead of in]) in this case we will get the expression 

J c 

We consider an isotropic environment in which the ten- _ -^i 

sor of conductivity is diagonal. In [3] the expression for the 
Laplace transform of the kind (fTl]) for the specific resis- 
tance p(w) is obtaned: Ml = ., . ■, {j;j)uj+ie;(c.oo) - d {J;J)u;c,t+{J;J)ij;c,t=l 

t[LL> + le) .V - ^ 



K 



P{^) - ^ = -^H^ + M]; (17) 

(7{L0) p[J,J) 

{J ; J)i^^ig 



Ki — — — — — ( j ; J)Lj+ie\(c,oo) ^ d ( j; J)uj-c,t- 
+ le) • / I 



{J,J) + {J;J)^^J ^''^ + (dc-f — 1 e-(J(c);J), 



i((jj + ze) 

If £ ^ 0, (r) ^ d ^ 0, and 

«^ ) J)u]-,{c,oo) 



where V is the volume of the system, (3 is inverse tem- 
perature, (7(w) is the scalar coefficient of conductivity. In 

the examples considered below in expressions (|17p - p^ it giwc^'j^^'j j-j _|_ ^j. J) ^ ^ 

is the value M that changes. Performing the operations of mi • • ^ n i r ■ i ^ tttti 

o j^l^jg expression at small values ol c is close to (llSp . I" or 

[3], with replacement of expression d by CH), in place of ^^^^ ^.^^ distribution p,{u) of the kind dn]) 

correlation (fT51) we obtain a more complicated expression ,1x1 j. c /Tttk • u ^-i, ^ j u ^.u 

^ — ' ^ ^ the Laplace transform p4|l is substituted by the opeation 

of the kind 

^^ / a \^ . . (^;^)^+*s;(o,c) + e^(^2-ei)(y^.5)^^^^^.^^_^^^ 

( J; J)^;a.fc + ^ (^-^ j ( J; J)^+.e;(c.oo) 

X ' where (A; B)^+.e;(o,c) = /dt e*(-+-)*(A(t), B), 



(A; i?}jj+,j£2:(c,oo) is given in and the value M from 

if=(J(0),J)--(J;J)^;„,fe+i- 

a psp IS replaced by 

(l- ^ e^"" f {J{c),J)+ ,v v^ . i*('^ + *ei),v ■r^ 

At e ^ and (r) ^ 00, A' ^ (J(0), J)--(J; J)^;a,fe+i+ = (J(0), J) - e'^^+'''^^{J{c); J) + ( j; J)^+.e,;(o,c)- 



a 

A; 

a 



(j;J).;a,fe + (^) (j;J).;(c.oo). gc(e.-e,)!0^i±^ [e^(-+-=)-(J(c), J) + ( j; J)^+,,,^(,,^) 

^^a + cJ 1{UJ + 182) L ' 
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At (r) oo and £2^0 the value M changes unessen- where k is wave vector, v{k) — 47re^/fc^, Si{k), the equi- 
tially, taking on the form hbrium statistical structure factor of ions, ee{k,Lu), equi- 



M = 



( j; j)uj+iei;{o,c) + 6 ^^1-^ — j; j)j^.((,_oo) librium dielectric constant of an electronic subsystem. 

If we pass from the correlation function (J; J)i^ of a 



Ks = (J(0), J) - e'("+'''i)'=(J(c), J) + ( j; J),^+iei:(o,c)+ kind (HH) to the correlation function of a kind (HH) in a 



e^-^(J(c),J) + (J;J}^;(e,oo) 



case 5aJ ^ at (F) ^ oo, e in ([8]) we see that at w = 0, 

11 1 1 ■ , ■ n 1 

where — is the expression tor the 



At small values of £i this expression is close to (fT8| . From a zub (a + c)*^ tr z-uf) 

conductivity , derived in [35] . For a case <?& ) with the 



(fT3l) it is seen that linie^^o £2(1") = e^^^''. ^ ^ 

distribution (u) (fTO|) at (F) — * 00, e — * 0, — = — . For 
Let us summarise explicit results for Coulomb systems. <J zub 

a case 3c) with the distribution at (/^} ^ 00,62 ^ 

Such systems were investigated by the NSO method in 

in (UnD, i = - e-^i"=. 
[MllSZllMllSlin]- We will follow [35]. We will derive the cr^"6 

Direct comparison of the obtained results to the ex- 
expressions for the conductivity of a completely ionised 

periment presents certain difficulties, since the parameters 

Coulomb plasmas in a constant electric field. For simplic- 

c,a,k,ei are apriori unknown. 

ity we limit ourselves to the case of the plasma consisting 
of two subsystems, electrons and positive ions. An isother- 
mal limit is considered when the characteristic thcrmaliza- 5 Superstatistics from distribution containing 

tion time for the charge carriers is much less than the re- 

^ lifetime 

laxation time of their composite momentum. The formula 

for the isothermal conductivity (at the frequency w = 0) In [JT] it is pointed out that a nonequilibrium distribution 
from [35] is written as follows: is characterized by an additional parameter related to the 

deviation of a system from the equilibrium (caused by the 

1 P / ^ \ ' ' 

— = --77 I I lim (J; J)ie, (19) field of gravity, electric field for dielectrics etc). In the 

where V is system volume, /3 = l/T is inverse temper- present work we consider open nonequihbrium, stationary 
ature, m = rrie and e are mass and charge of electron, systems, certain point of metastable states. Investigations 
n = ne is average electron density. For the correlation on spin glasses and other aging systems, where a "waiting 
function from (HU) in an expression is obtained [35]: time" plays an important role, allows to anticipate the us- 

ability of this approach with respect to them as well. In 

. . the present paper we suggest a new choice of an additional 
lim (J; J)i^ = 

2 rill parameter in the form of the lifetime of a physical system 

--(-) Skk-'v{k)S,{k) hm-Im— — , 

P ^fn/ \u>^o Lu ef.{K,,u!)\ which is defined as a first-passage time till the random 



10 V.V. Ryazanov: Lifetime distributions in the methods of non-equihbrium statistical operator and superstatistics 

process y{t) describing the behaviour of the macroscopic The distribution ()20p with the hfetinie contains two 

parameter of a system (energy, for example) reaches its different time scales: the first relates to the energy E, 

zero value. The lifetime (or F) is thus a random pro- and the second - to the lifetime itself F, this latter one 

cess which is subordinate (in terms of the definitions of accounts for large-scale time correlations and large-time 

the theory of random processes [42]) with respect to the changes in E by means of a thermodynamic conjugate to 

master process y{t), the lifetime value 7. The similar operation can be derived 

starting from NSO. The structure factor uj{E) is thus re- 

F,^znf{t:y{t)^0}, y{0) ^ x > . 

placed by lo{E, F) - the volume of the hyperspace contain- 

This definition of the lifetime is taken from the apparatus -^^ ^-^^^ ^^^^^^ of Es^ndF. The number of phase points 

of the theory of random processes where it is widely used between {E, E + dE; F,F + dF} equals cj{E, F)dE dF. 

in the theory of queues, stochastic theory of storage [SI, The value thermodynamically conjugated to the lifetime 

Kramers problem of the escape rate out of a potential well ^^^^^^^ ^he entropy fluxes and entropy production 

nags] and so on. These questions are discussed in text- ^^ich characterize the peculiarities of the nonequilibrium 

books by van Kampen gS], Gardiner 07] and many other processes in an open thermodynamic system. If 7 = 

[IH]. The lifetime plays part in the theory of phase tran- ^^^^ /3 = /3o - l/ikBT^q), where ke is the Boltzmann 

sitions, chemical reactions, in the dynamics of complex constant, T^, is the equilibrium temperature, then the ex- 

biomolecules etc. pressions (UHll-dll yield the equilibrium Gibbs distribu- 

Using a maximum-entropy principle gH], it is possible ^^^^^^ Que can thus consider (I2ni)-(E21) as a generalization 

to derive the form of the expression for microscopic (but ^j^g Gj^bs statistics towards the nonequihbrium situa- 

coarse-grained) probability density in the extended phase ^j^n. Such physical phenomena as the metastability, phase 

space [TSl[Tnj[2ni[lI] transitions, stationary nonequilibrium states are known to 

, „ r ir7/n \ violatc thc eQuiprobabilitv of thc phasc spacc poluts. Thc 

p{z;E,F) ^exp{-PE --fF}/Z{f3,-f), (20) 

value 7 can be regarded as a measure of the deviation 

where 

from the equiprobability hypothesis. In general one might 
Zi0.,)^leM-^E-,n4. ^ (21) choose .h. value r .s . subp„ce. of .o™ o.he, W 

as chosen above. Mathematically the introduction of the 

lifetime means acquiring additional information regarding 

is the partition function, 13 and 7 are Lagrange multipliers j 1 • j_ i j_- 1 1 • 

^ If/ " " an underlying stochastic process; namely, exploring the 

satisfying the equations for the averages / ^ i- \ ^- r ^ u j 1 

J o ^ t. (stationary) properties of its slave process beyond merely 

/z7\ d\nZ dlnZ knowledge of its stationary distribution. 



dE dFuj{E, F) exp{-PE - -fF} 
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In the distribution ((20|) containing Ufetime as a ther- and so on. 

modynamic parameter, the probability for E and r is We note that there is a similarity between the method 

equal to of superstatistics, where the averaging is performed over 

Pi^T^) = ■ (23) a parameter /? (for example, the inverse temperature), as 

TT • • , ,1 j77?ni n 1 i • ii 1- i -1 in (HI), and the method of NSO, where averaging is per- 

Havmg mtegrated (|23|) on 7 , we obtam the distribution " 

„ , . , formed over the extension of past time u — t ~ tn, as 

of a kind 

^ in J!]) [51]. Expressions H]), ^ and ^ are described 

P(E)= I P(E,r)dr = [ uj{E,r)e-^^dr. , ^. , , , , rr^ rr., , 

J Z((3,'j) by the subordinated random processes j42j. Ihe Zubarev 

(24) 

approach claims that the source term should be infinites- 

The structural factor uj(E,r) has a sense of the joint . „ n mi • i ;i • i ■ 

imally small, ihe question is now whether a vanishing 

probability for £' and considered as a stationary distri- ^ ii • ii i-rc 

source term would yield results dinerent trom the super- 

bution of this process. We shall write down ^ ^. ^. . i r 

statisrtics as well. In [1611171152] a very simple example of 

n 

t -n\ t ^ r^^ ^ Ip^ d y ^ ir\ ^ot;^ the Browniau particle is Considered. Its velocity u Satisfies 
u(E,r) ^ uj[E)uJi{E,r) ^ uj{E) 2_^Rkfk(r,E) . (25) ^ 

k=l 



In the last equality ([^5)) it is supposed, that there exist n , , 

is Gaussian white noise, 7 > is friction constant, and 



the linear Langevin equation v — —71; + aL{t) where L{t) 

classes of ergodic states in a system; Rk is the probabil- , , r , • • ,1 , • 1 r , • 

the strength of the noise is controlled m a usual fashion 

ity that the system is in the fc-th class of ergodic states, , , . , , •,■ 1 • r 

by the parameter a. ihe stationary probability density of 

fk{r,E) is the density of lifetime distribution F in this . „ . . , , . „„i , 

V is Gaussian with average and variance p , where the 

parameter /3 = can be identified with the inverse 



class of ergodic states (generally fk depends on i?). As a 2 



physical example for such situation (typical for metals or ...... 

temperature of the statistical mechanics (we assume that 

glasses) one can mention the potential of many complex 

the Brownian particle has a unit mass). This simple situ- 

systems. This case is considered in |50|. . . , • , 

ation completely changes if the parameters 7 and a m the 

In [2U1[^ various models of superstatistics are ob- , . . , . , n 

stochastic ditterential equation are assumed to fluctuate 

tained from ([20]) -(gS]). For example: „ , -n w j ^ ^1 . -^1 

as well, io be specinc, let us adopt that either 7 or cr or 



op 

P(^jt;^ ^ ^ . both fluctuate in such a way that /3 — "f/a^ is x^- dis- 



1 + (g - 1) ( 1 - e-PoroE 



tributed with degree n. This implies that the probability 
density of (3 is given by 
P{E) - „ 

^^^> r(n/2) {2P0J ' ^\ Wo 



q-l 



l + (g-l)(l- [l + (gi-l)(l-e-/5oroi?j]«i-i)l (5^ ^ j Pf{P)dP- (26) 
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In the Zubarev approach if the parameter e — > in the 
exponential distribution (|4]), the source in the LiouviUe 
equation vanishes. Relating to the distribution (|26p it cor- 
responds to /3o ^ oo and (n/Po) — > 0. As w <-> /3, Pq{u) ^ 
/(/3), (r) ^ /3o in HI), Q, (HSl), the case (F) oo cor- 
responds to that /3o — > oo. Apparently, this is the limiting 
case of CT'^ — » 0, when no stochastic element is present, 
but the system is subject to the dynamical force only, as 
in the LiouviUe equation without random source. 

6 Piecewise continuous distributions for 
functions R, f from correlations (l^4 D -([^5l ) . 

m 



fiix,E) 



o-ir 



f{x,E,r)dr. 



We note, that the correlation (P7)l includes Laplace trans- 
form to which a probabilistic sense can be ascribed accord- 
ing to ^ . For / and /i from ([27]) it is possible to use the 
models P^[T^[2U1I^ which leads to superstatistics of the 
kind exp{— yi?}. The similar approach developed in [7j al- 
lows to obtain the correlation value /o(j/) ~ ae^^^ for the 
model of phase synchronization, which under certain con- 
ditions reduces to the form exp{— yi?}. The parameters a 
and k depend on the problem. 

Let us consider this problem for a simple case of the 
function R{y) from (|27p . allowing to write an obvious form 
of probability density where the function R{y) represents 



a combination of delta-function and homogeneous distri- 
In this section the distributions of lifetime, having a differ- bution: 



Riy) = <^ 



Then 



ent shape on the different temporal intervals of evolution 
of the system, are considered. Such behaviour is charac- 
teristic for many physical systems. It is stressed in [29) 
that non-equilibrium systems can have different stages of 
evolution. In 53J it is shown that the first passage time 
probability density distribution changes depending on the 
value of control parameter. Solutions of the Kramers equa- p[E) — j R{y) ^ 
tion, which are related to the first passage time probability 
density distribution also depend on the control param- 
eter. Such transitions in real systems are widely encoun- 
tered, the aging of materials is just one of known examples 
thereto. 

Let us write the expressions in the form 

oo 



pS{y-a), y<c, a<c, p<l; 
(1 — p)m^^, 0<c<y<c + 'm; 
y > m . 



dy 



P 



(1 + -fae'^y^) 



( 

+ (l-p) 



1 

'kE 



■In 



V 



okEc 



7a 



7a 



Let us choose now the expressions of a kind 

P{E). 



R{y)\e~''y^ ^y 



z 



P{E) = 



R{x)h{x,E)dx- 



(27) 



(coinciding with ([T])) with piecewise continuous distribu- 
tion of the function R{y). If we set the function R{y) in 
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the form of the gamma-distribution with different values 
of the parameters a and r in different areas, 
A" 



R{y) 



5 ^ „,r-lo-A£ 



r(r) 



y 



y <c; 

y>c, 



, , rir)r{a,\c) 

then b = „ , { „ , . and 



piE) = 



r{a)r{r, Ac) 
1 



r{a, Ac) 



7(a, (k^; + A)c) 
r{r, {kE + A)c) 



(1 + kE/X)"' 
1 



(28) 



r(r,Ac)^ ' "'^'{1 + kE/XY 

This distribution is more complex, than Tsallis distribu- 
tion [14ill5j . obtained from gamma-distribution by means 
of a method of superstatistics [TBlfTT] . Multipliers in a form 
1/(1 + kE/X)" correspond to the Tsallis distribution, but 
in other factors depending on E are present. The dis- 
tribution (pS)) passes in 1/(1 + kE/X)°' at r ^ a. 
If 

A" 



R{y) 



r{a) 



I (y + a) 



k+l ■ 



y > C, 



thcn6^^(" + -)^^("'^-) and 



p{E)^ 



r{a) 
1 



zr{a) 



-f{a,{KE + X)c) 



{I^^kE/XY 



k{a + c)*r(a, Xc)e'"'''{KEYr{-k, kE{c + a)) 

Other similar examples can be considered as well. One 
could set three and more areas of piecewise change of vari- 
ables (for example, two areas corresponding to different 
phases and a transitive layer between them). A contin- 



function or to be set by functions of a kind R{g{-)). A com- 
bination of discrete Rk and continuous R{y) distributions 
in different areas is also possible. 

Let us consider two more simple examples of the com- 
bination of gamma-distribution with delta-distribution and 
with homogeneous distribution: 

R{y) = I r{aY 



bS{y — d), y>c,d>c. 



Then h ■ 



p{E) 



r{a, Ac) 
1 



zr{a) 



7(a, (kE + A)c) 



(1 + kE/XY 

r{a,Xc)e-^^^ 



By 



r A« 



R{y) = < 



r(a)- 

bra^^ , 0<c<j/<c-|-m; 
0, y > m , 



r{a, Ac) 



p{E) = 



zr{a) 



7(a, [kE + X)c) 



1 



{1 + kE/XY 



r{a,Xc)e-^^^{nE)-^{l 



Various combinations of functions of distribution from 
uous change of parameters can be treated on equal foot- |161I17| can be brought into consideration: for example, 
ing, for example, considering the parameter of the gamma- lognormal superstatistics, gamma superstatistics and in- 
distribution continuously changing with some distribution verse gamma superstatistics. 
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7 Applications 

Passing to superstatistics we get different distributions 
R{y) from ^7} at y < c and y > c or f{l3) from ^ 
at P < c and (3 > c. For example, exponential or gamma- 
distribution at /3 < c and Pareto distribution correspond- 
ing to Tsallis distribution at /3 > c. A multitude of vari- 
ous combinations of different distributions can be encoun- 
tered, including, for example, gamma-distribution with 
various parameters at different temperatures, the subordi- 
nated distributions and so on. The parameters of gamma- 
distribution can change discretely, but can be continuous 
as well. 

In [5S] it is assumed that Beck and Cohen's super- 
statistics provides a suitable description for systems with 
mixed regular-chaotic dynamics. Such systems can be de- 
scribed by means of the suggested approach. Examples of 
different behaviour of systems at different temperatures 
are obvious including superconductivity, and superfluid- 
ity, and other phase transitions. 

The approach which is more general, than the super- 
statistics theory, consists in setting piecewise continuous 
distributions for R{y) from expression (|27p . that is the 
probabilities for a system to be in in the fc-th state. 

Such distributions can describe laminar and turbulent 
modes in the stream, for example, the diffusion of the to- 
bacco smoke flow in the atmosphere. For this case the 
index k (the parameter ?/ - in a continuous case) corre- 
sponds to the spatial coordinate of the flow. In any point 
c there is a transition from a laminar mode to turbulent. 
The distribution R{y) or f{(3) can be described by the 
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correlations obtained in '56] (lognormal distribution) . The 
situation described here is more general, than the only one 
transition at a certain temperature, as in the case of su- 
perstatistics. 

In [ST the entropic index can take various values, as 
in the present work. In [58j the entropic parameter q of 
the Tsallis distribution depend on the parameter m, and 
the parameter m is used to give account of an energy loss 
rate or energy dissipation rate (or perhaps, the energy ab- 
sorption rate). The Tsallis distribution appears to depend 
on the parameter m. Such situation corresponds to the 
approach of the present work. 

The distribution suggested in [59| represents a spe- 
cial case of the piecewise continuous distribution used in 
the present work. In [BD] the problems close to the prob- 
lems of the present work are considered: the occupation of 
the accessible phase space (or of a symmetry-determined 
nonzero- measure part of it), which in turn appears to de- 
termine the entropic form to be used. 

The results similar to results of Section 5 are obtained 
in [61| . For example, the distribution /(pi . . . pn) of 
metastable states with local fluid densities pi in different 
spatial domains i — I . . . N depending from the exponent 
A is again related to the distance from the conventional 
equilibrium, as does the value 7 from Section 5. 

Besides the spatial heterogeneity the piecewise con- 
tinuous distribution can describe the time changes. The 
suggested approach allows to use the methods of the the- 
ory of random processes for treating specific problems; for 
example, to refer to the stochastic theory of storage [131 
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I62j . setting rates of an input in a potential well and an 
output from it, and to the methods of the Kramers escape 
problem [iilUS] . 

8 Conclusion 

As it is stated in [53] , the existence of different time scales 
and the flow of the information from slow to fast degrees 
of freedom create the irreversibility of the macroscopical 
description. The information thus is not lost, but passes in 
the form inaccessible at the Markovian level of description. 
For example, for the rarefied gas the information is trans- 
ferred from one-particle observable to multipartial correla- 
tions. In [4J the values e = 1/ (r) and Pq{u) — eexp{— eu} 
are expressed through the operator of entropy production 
and, according to the results of 63J, in terms of the flow of 
the information from relevant to irrelevant degrees of free- 
dom. The introduction of the function Pq{u) in NSO cor- 
responds to the specification of the description by means 
of the effective account of communication with irrelevant 
degrees of freedom. In the present work it is shown, how 
it is possible to expand the specification of the description 
of memory effects within the limits of NSO method. A 
more detailed description of the influence of fast varying 
variables on the evolution of system is suggested based 
on specifying the density of the life span distribution of a 
system. 

In many physical problems the finiteness of the life- 
time can be neglected. Then e ^ 1/ (F) 0. For example, 
for a case of the evaporation of liquid drops it is possible 
to show 64], that non-equilibrium characteristics depend 
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on exp{y^}; y — e/(2A2)^/^, A2 is the second moment 
of the correlation function of fluxes averaged over quasi- 
equilibrium distribution. Estimations show, that even at 
the minimal values of lifetime of drops (generally of fi- 
nite size) the maximum sizes are y — 6/(2X2)^^^ < 10^^. 
Therefore finiteness of values of (F) and e does not influ- 
ence the behaviour of system and it is possible to consider 
£ = 0. However in some situations it is necessary to take 
into account the flniteness of lifetime (F) and values £ > 0. 
For example, this is the case of the nano-drops. 

Changes of the form of the source in the Liouville equa- 
tion, as well as the expressions for the kinetic coefficients, 
average fluxes, and kinetic equations can be obtained with 
the use of the NSO. It is possible to choose a class of life- 
time distributions for which after thermodynamic limiting 
transition and tending the average lifetime of system to in- 
finity the results are reduced to those obtained under ex- 
ponential distribution for lifetime, used by Zubarev. How- 
ever there is also another extensive class of realistic distri- 
butions of lifetime of system for which even if the average 
lifetime of system tends to infinity the non-equilibrium 
properties essentially change. It is a consequence of the 
interaction of the system with its environment. 

For the distributions of the kind (dH]), (HH), having 
different form for different argument spans, non- vanishing 
corrections to Zubarev NSO persist even for infinitely 
large systems with infinitely large lifetimes. In the present 
work it is shown, that this situation is possible, for exam- 
ple, for the distributions of lifetime of a system, having 
different form at different stages of the evolution of a sys- 



It is interesting to investigate possibilities of such a 
choice of the function pq (u) which would in most full fash- 
ion correspond to physical conditions in which a system is 
placed. For an optimal choice of a form of function Pq{u) 
it is possible to use the method of maximum entropy prin- 
ciple. 
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tem. Such behaviour corresponds to the realization of the statistics when /3o ^ oo, and the distribution of an inten- 
evolution of a system in a number of subsequent stages sive parameter has a piecewise continuous character. 
|29| . The form of distribution chosen by Zubarev for the life 

span represents a certain limiting case. The choice of the 
lifetime distribution in NSO is related to the account of 
the past of a system, its physical features, on the present 
moment; for example, with the account of the age of a 
system only, as in Zubarev form of NSO [T1[51[31|311S1[3] at 
e > 0, or with more detailed characteristic of the past 
evolution of a system. The obtained results are essential 
The use of the theory of the superstatistics in var- in cases when it is impossible to neglect the memory effects 
ious applications is related to the piecewise continuous since the memory correlation time there is not vanishing, 
functions for the density of distribution of lifetime of the The analysis of the corresponding time scales is necessary 
system. The task of these functions of distribution allows as it is noted in |63j . 

to receive new expressions for the superstatistics, corre- The main objective of the present article is to show, 
sponding to various physical situations. Similar results how the systems with infinitely large average lifetime can 
seem to be useful, for example, in the investigation of small induce nonvanishing sources in the Liouville equation, and 
systems. A number of results following from the interprc- in what consequences for the method of the NSO it results, 
tation of NSO and Pqiu) as a density of lifetime distri- Superstatistics with piecewise continuous distributions of 
bution of system [5], can be obtained from the stochastic intensive parameter are considered as well, 
theory of storage [43J and theories of queues. For exam- 
ple, in 031 the general result that the random variable of _ _ 
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